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^ ■ Abstract 

h^ , We perform the phase space analysis in terms of the linearization technique in the Hoi^ava-Lifshitz 

"^ I gravity with the softly broken detailed balance condition. It can be shown that the bouncing 

universe appears only for the positive spatial curvature of A; = +1, and it is possible to obtain 
oscillating universe with the help of the negative dark radiation and the negative cosmological 
constant. 
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Recently, Hofava-Lifshitz (HL) gravity motivated by the Lifshitz theory in the condensed 

matterphysics [l| has been proposed as an ultraviolet (UV) completion of general relativ- 

^. The basic ingredient of the UV completion comes from an anisotropic scaling 



ity 



between space and time, 



t^bH, 



x* -)■ 6x* 



(1) 



where the conventional form can be recovered for the infrared (IR) limit of 2; = 1. In the 
original HL gravity, the "detailed balance" condition (DEC) was introduced to simplify its 
quantum behavior and renormalization properties, and it has been shown that its renormal- 
izability depends on the renormalizability of the topologically massive gravity (TMG) 5|, |6| 
used for the DBG [7|, |8|. However, it is well-known that HL theory with the DBG can only 
admit solutions which are asymptotically AdS. To solve this issue, Hofava considered a softly 
broken model which contains IR modification without any change in UV terms, so that HL 
theory admits asymptotically flat solutions as well without loss of the renormalizability 3|. 

As pointed out by Hofava in Ref . |3| , the softly broken HL theory still satisfies the DBG 
at short distances, since the dominant terms (sixth order terms in spatial derivatives) are 
constructed by the Gotton-York tensor through the DBG. The Gotton-York tensor can be 
derived by a variational principle with the gravitational Ghern-Simons (GGS) action p, |6|. 
In other words, the renormalizability of the softly broken HL theory in the UV limit, depends 
on the renormalizability of the GGS term. Even though the subdominant terms (fourth order 
terms in spatial derivatives) are taken into account, they can be constructed from the TMG 
with the vanishing cosmological constant, so that its renormalizability depends on that of the 
TMG. Therefore, it means that the present softly broken HL theory is still renormalizable. 

On the other hand, there has been much attention to HL gravity to explore the cosmo- 



logical implication of HL theory |9l-l25| and related subject 26l-l3l|. As an extension of HL 
theory, F{R) HL gravity has been well appreciated to show that an accelerating phase of 



the universe can be obtained 



32l |33| , and a unified description of inflation with dark energy 



is possible 3^, |35|- In particular, bouncing universe has been intensively studied 36|-|45| in 
general relativity (for a recent review, see Refs. |46l. l47l| ) . 

In this note, we would like to perform phase space analyses of the HL gravity with the 
softly broken DBG. In fact, the original DBG severely restricts the asymptotic behaviors 
so that the only asymptotically anti-de Sitter spacetime is possible. What it means is that 
at the late-time stage of the universe the decelerated expansion may appear, unfortunately. 



which seems to be incompatible with the recent accelerated expansion of the universe. By 
the way, it is possible to obtain the asymptotically flat solution of the vanishing cosmological 
constant along with the asymptotically anti-de Sitter solution for the softly broken DBG. So, 
one may ask what the crucial cosmological difference is between the case of the vanishing 
cosmological constant and the nonvanishing one in the softly broken DBG model, which 
is a main motivation of the present work. Eventually, it turns out that the cosmological 
behaviors between them are not so much different from each other. It shows that only for 
k = +1, cyclic universes exist throughout matter bounce at early stage and decelerated 
expansion at late time for both cases. As for k = —1, there exists an unstable fixed point in 
the forbidden region of the negative total energy for both cases. In this case, we can avoid 
the initial singularity by imposing the positivity of the total energy. Unfortunately, all of 
them do not show the recent accelerated expansion, which is a weakness of the model with 
the DBG either with the softly broken DBG. 

Now, considering Arnowitt-Deser-Misner (ADM) decomposition of the metric with ds'^ = 



-N'^c^dt'^+gij{dx' + N'dt){dx^ + N^dt) |48|,|49|, the Einstein-Hilbert action can be rewritten 
as 

Ieh = j7^ I dW^ [K - 2A] 
IbTrGAf J 



/ dtd^x^N [KijK'^ -K^ + c^{R~ 2 A) 



IGvrGjv 
where Kij = -^^ [gij — ViNj — V jNi] is the extrinsic curvature at t = constant hyper- surface, 

and the dot denotes the derivative with respect to time t. Here, gij, R, and Vj are the metric, 

the intrinsic curvature, and the covariant derivative in the three-dimensional hyper-surface, 

respectively. The above anisotropic scaling ([T]) gives A^* — > 6^~^A^* and c — t- b^^^c while gij 

and A^ are invariant. Requiring that the Planck constant is also invariant, h -^ h, we get 

E —)■ b^^E and M -^ b^^^M, where E and M are energy and mass, respectively. Then, the 

Hofava-Lifshitz gravity (HL) under the softly broken DBG [3] is given by 

Jhl = fdtd'x^Nl^ {K,,K^^ - XK^) - $-^C,fi'^ + '^e'^^ R,,V ,Ri 

where k"^ is a coupling related to the Newton constant Gn, and A is an additional dimen- 
sionless coupling constant. In addition, u represents the IR modification which is essential 
to have asymptotically fiat solutions. The couphng constants fj,, Aw, and ( come from the 



m 



three-dimensional Euclidean topologically massive gravity action [5|, l6| given by 

1 



W 



/i / d'x^{R - 2Ah^) + 



x(r), 



(4) 



where x(r) represents the gravitational Chern-Simons term. Then, the scaling of couplings 
can be obtained as k^ — )■ b^'^n"^, jj, -^ b~^fi, Aw — )■ b~'^Aw, C ~^ C and u — )■ b~'^u. Note that 
identifying the fundamental constants with 

^2 /;x2(^_A^^) 



K 
T 



G 



K^C^ 



Af 



A 



3A2 



VK 



(5) 



3A-1 ' " 327r' " 2(a;-Aw', 

the Einstein-Hilbert action ([2]) can be recovered in the IR limit when A = 1. For A > 1/3, 
we have u > Aw and A < from the identification ([5]), while for A < 1/3, we have u < Aw 
and A > 0. From now on, we will assume A > 1/3 to keep the Einstein limit. 
Now, we are going to consider the Robertson- Walker (RW) metric, 

dr'^ 



ds^ = -c^dV + a\t) 



1 — kr'^ 



+ r'^dnl 



(6) 



where fc = 0, ±1 are normalized spatial curvature. After some tedious calculations, the 
equations of motion from ([3]) can be obtained as 



3(3A-l)/72 



{?>\~l)[H+-H^ 






A 
3" 



k 



k"^ 



2a^{uj-Aw)_ 



Y [P + Pvac + Pk + Pdr] , 
A k 

T"^ 3^ 



K 



— V — 3c 



e 



6a^{io - Aw) 



K 



- — [p + Pvac +Pk+ Pdr] , 



(7) 



(8) 



where H = a/ a is the Hubble parameter. Apart from the normal matter contributions of 
p and p, the additional energy-momentum contributions from the potential terms ([3]) are 
explicitly written as 



12c2A 



Pv 



-Pv 



Pk 



1 

-gP. 



4c2 k 



9 9 ) Pdr 

K^ a^ 



1 

gPdr 



4c2 



fc2 



(9) 



fi;2 1a\uj - Aw] 

Similarly to the general relativity, p^^c and p^ come from the vacuum energy and the spatial 
curvature contributions, respectively. In particular, in the HL theory, p^r is called the dark 



radiation. The total actual energy can be defined by ptot = P • 
negative dark radiation is enough to realize the matter bounce 



Pva 



13 



Pdr- Note that the 



14| . so we do not need 



any additional ad hoc energy sources to make the bouncing universe. By the way, combining 
Eqs. ([7]) and ([8]), we can get the acceleration of the scale factor, 



a K^ , „ , K^ 



(3A - 1)- = -- [ptot + Sptot] = -Y^[l + 3c^(«)]ptot, (10) 

where we used the relation p^ + 3pk = and u^a) = Ptot/Ptot is the equation-of-state 
parameter for the total energy. 

Now, let us perform numerical analysis for this system ( ITOl) and obtain phase portraits. 
First of all, we should introduce additional variable f = a, so that we get two first-order 
differential equations, 

a = V, (11) 

.^-2-Jl,Ma)]U^^), (12) 



where we used Eqs. ([7]) and ( ITOj) . Next, following the linearization technique [50|, the 
Jacobian matrix is obtained as 



A 



1 

-3c^kuj'{a*)/{3\-l)a* 



(13) 



and the corresponding eigenvalues are t] = ±^^/ —Sc^kuj' (a*) / (3X — l)a*, then fixed points 
(a*,0) are classified as a center for ku'^a*) > or a saddle node for ku\a*) < 0, where a* 
satisfies 1 + 3a;(a*) =0. 

We are now in a position to specify the matter source which consists of the conventional 
cold matter and radiation, p = Pm + Pr with pm ~ a~^ and pr ~ a^"^. Then, the equation- 
of-state parameter a; (a) is given by 

3 [ridr + fir + a^m + «^fivac] ' 

where the density parameters are defined as fij = Pi/pc evaluated at the present universe 
scale ao, which can be fixed to ao = 1 for convenience. Assigning fivac — —0.7, Qm — 0.3, 
Qk ~ — 10~^/c, Qr ~ 10~^, and Q^r ~ —1.2 x 10~^, the phase portraits for this system are 
plotted in Fig. [T]in the logarithmic scale for (a) k = +1 and (b) k = —1, respectively. The 
black and white dots on the a-axes represent fixed points. It is of interest to show that in 
Fig. [T]the bouncing universe appears only for the positive spatial curvature oi k = +1, and 
it is possible to obtain oscillating universe with the help of the negative dark radiation and 



(a)k=+l 



(b)k=-l 





ln(a) 







ln(a) 



FIG. 1: The phase portraits are plotted in the logarithmic scale for (a) k = +1 with the stable 
fixed point and (b) k = —\ with the unstable fixed point. The black and white dots on the a-axes 
represent fixed points (a*, 0) satisfying 1 + 3a;(a*) = 0. 



the negative cosmological constant. We have maintained the positivity of the total energy 
density. In this respect, the inside between the horizontal arrows in Fig. [1] (b) is not allowed. 
As for the negative contribution of dark radiation, it is interesting to note that there have 
been concrete justifications for models with negative density, in particular, a brane universe 
moving in a curved higher dimensional bulk space [51| and a model of dark energy stemming 
from a fermionic condensate 52|, although we do not know how to justify the negative dark 
radiation in this model. In Fig. [2l the cosmological evolution can be shown for the case of 
the vanishing cosmological constant following the same procedure. Eventually, it turns out 
that the cosmological behaviors between the cases of the vanishing cosmological constant 
(Fig. [2]) and the nonvanishing cosmological constant (Fig. [1]) are not so much different from 
each other. 

Final comment is in order. Unfortunately, the HL cosmology with the DBC does not 
exhibit the recent accelerated expansion of our universe. On the other hand, a positive 
cosmological constant can be obtained for A > 1/3 by an analytic continuation, /i — )■ i/i and 
(^^ — > —iC,"^ J26|, |27|, which makes the coefficients of the potential terms in the action ([3]) 
have the opposite sign, so that we have u < Kw and A > 0. Then, it gives the late- 
time accelerated expansion, however, it causes the initial singularity problem because of the 



(a)k=+l 




ln(a) 




/ 







^y ln(a) 



FIG. 2: The phase portraits with same data as in Fig. [T] except for vanishing cosmological constant 
are plotted in the logarithmic scale for (a) k = +1 with the stable fixed point and (b) k = —1 with 
the unstable fixed point. The black and white dots on the a-axes represent fixed points (a*,0) 
satisfying 1 + 3oj{a*) = 0. 

absence of the matter bounce. 
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